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Abstract 



In the framework of Euclidean QCD on a torus, we study the spectrum 
of the Dirac operator through inverse moments of its eigenvalues, averaged 
over topological sets of gluonic configurations. The large- volume dependence 
of these sums is related to chiral order parameters. We sketch how these 
results may be applied to lattice simulations in order to investigate the chiral 
phase transitions occurring when Nf increases. In particular, we demonstrate 
how Dirac inverse moments at different volumes could be compared to detect 
in a clean way the phase transition triggered by the suppression of the quark 
condensate and by the enhancement of the Zweig-rule violation in the vacuum 
channel. 
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I. INTRODUCTION 



Understanding the Spontaneous Breakdown of Chiral Symmetry (SB^S) remains one of 
the most challenging non-perturbative problems of QCD. Forthcoming experiments [l] [| 
should reveal some of its features, at least in the non-strange sector in which the effective 
number of light quark flavours is minimal (Nf = 2). It is generally expected that if Nf 
increases (keeping the number of colours N c fixed), the theory meets phase transitions and 
the chiral symmetry is eventually restored. The argument is originally based on properties 
of the QCD /3-function in perturbation theory. The well-known statement of the "end 
of asymptotic freedom" for Nf > WN c /2 || has been further completed by the analysis 
of the so-called "conformal window" || suggesting a restoration of chiral symmetry for 
lower Nf, such as Nf ~ 10 (for N c = 3) ||. Less perturbative and more model-dependent 
investigations, based on a gap equation 0] or on a "liquid instanton model" || , also indicate 
that a chiral phase transition could occur for Nf substantially below lliV c /2. 

It is important to understand, at least qualitatively, the non-perturbative origin of the 
suppression of chiral order parameters for an increasing Nf. We have recently argued that 
such a suppression might result from a paramagnetic effect of light (massless) quark loops 
[T0|1 , i.e. it could be due to "sea quarks" and, consequently, it could escape a detection in 



quenched lattice simulations, or in any other approach neglecting the fermion determinant. 
An appropriate framework to develop these ideas and to ask precise questions is the formu- 
lation of QCD in an Euclidean box L x L x L x L, with periodic (antiperiodic) boundary 
conditions for gluon (fermion) fields, up to a gauge transformation. In this framework, the 
SB^S pattern is reflected by the dynamics of lowest eigenvalues of the Dirac operator: 

H[G]= 1 ,(d, + iG,). (1.1) 

This Hermit ian operator has a symmetric spectrum with respect to zero: {if, 75} = 0. 
Positive eigenvalues X n are labeled in ascending order by a positive integer (one further 
denotes A_ n = — X n and 0_ n = 75 <fi n for the corresponding eigenvectors). SB^S is related to 
a particularly dense accumulation of eigenvalues around zero ]TT|- p^ , pT| ] . Models of such an 
accumulation in terms of random matrices [14] or instantons [|U|] have been proposed. Some 



chiral order parameters are entirely dominated by the infrared extremity of the spectrum 



of the Dirac operator ( |1 . 1| ) . This makes them particularly sensitive to the statistical weight 
given to smallest Dirac eigenvalues in the functional integral, which is suppressed in the 
massless limit by the iVj-th power of the fermion determinant. A good example is the quark 
condensate, defined by: 

£(#/) = - lim <0|uu|0), (1.2) 

mi,m2,—mtf,— >0 

where 7%. . . mjv, denote the Nf lightest quark masses and u represents the lightest quark 
field. H(Nf) receives exclusive contributions from the smallest Dirac eigenvalues that behave 
in average as 1/L 4 , and it is consequently expected to be the most sensitive order parameter 
to the variation of Nf and to a phase transition. Other order parameters are less sensitive, 
like F 2 (Nf), defined as the SXJi(Nf) x SUr(./V/) limit of the coupling of the Goldstone bosons 
to the axial current: 
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F 2 (N f )= lim F*. (1.3) 

F 2 (Nf) may be non-zero due to Dirac eigenvalues accumulating as 1/L 2 ||13|| . For this 
reason, F 2 (Nf) should exhibit a weaker ^/-dependence than E(iV/). Finally, observables 
with no particular sensitivity to the infrared edge of the Dirac spectrum (p-mass, string 
tension, etc. . . ) have no reason to be strongly affected by the fermion determinant and by 
the iVj-dependence. 

Let us first consider the thermodynamical limit and denote by n cr i t (N c ) the critical value 
of Nf at which the first chiral phase transition takes place. Just below n cr i t (N c ), the order 
parameter £(iVy) drops out, whereas its fluctuations may be expected to become important. 
We have shown || that the latter would manifest itself by an enhancement of the Zweig- 
rule violation just in the vacuum channel J PC = ++ . An important Zweig-rule violation 
is precisely observed in the scalar channel |16|], and nowhere else (with the exception of 



the pseudoscalar channel driven by the axial anomaly). Whilst the signature of a nearby 
phase transition is rather clear just below n CI it(N c ), it is more speculative and ambiguous 
above the critical point. First, above n CTit (N c ), colour might still be confined (confinement 
has no obvious relation to small Dirac eigenvalues). Second, despite T,(Nf) = 0, the chiral 
symmetry need not be completely restored. The Goldstone bosons coupling to conserved 
axial currents with the strength F(Nf) might survive to the iV/-induced phase transition. 
This is reflected by the possibility that the iV/-sensitivity and suppression of the order 
parameter F 2 (Nf) might be considerably weaker than in the case of the quark condensate 
|T3|j . Of course, this is a highly non-trivial possibility, which presumably depends on the 



existence of a non-perturbative fixed point in the renormalization group flowQ. Here, we take 
as a working hypothesis that above n CT i t (N c ), a partial SB^S still occurs, due to F 2 (Nf) ^ 0. 
The results of our paper allow, in particular, to test this hypothesis. 

The central question remains how far is n cr i t (N c ) (for N c = 3) from the real world, in 
which the number of light quarks hardly exceeds Nf = 2 — 3. Some recent investigations 
actually indicate that n cr i t (3) could be rather small, and/or that the real world could already 
feel the influence of a nearby phase transition. First, some lattice simulations with dynamical 
fermions observe a strong iV/-dependence of SB^S signals for Nf as low as 4-6 |18| , [19| . 
Second, a method based on a well-convergent chiral sum rule has been proposed, which allows 
to study phenomenologically the variation of H(Nf) for small Nf | 20| . It has been found that 



existing experimental information on the Zweig-rule violation in the scalar channel leads to 
a large reduction of T,(Nf) already between Nf = 2 and Nf = 3. 

The purpose of this paper is to analyze in a model-independent way how iV/-induced 
chiral phase transitions manifest themselves in the finite-volume partition function. In par- 
ticular, we shall investigate the volume dependence of the inverse spectral moments of the 
Dirac operator (|1.1|): 



would imply 



1 If one sticks to cut-off-dependent bare quantities, it is possible to argue that X 
F = 0, i.e. the full symmetry restoration [17]. This argument is however based on an inequality 
for which it is by no means obvious that it survives in the full renormalized theory. 
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(L4) 

averaged over topological sets of gluonic configurations. For Nf <C n crit (N c ), the leading 
large- volume behaviour of such inverse moments has been worked out in details by Leutwyler 
and Smilga |21] . In order to investigate how this result is modified in the vicinity and above 
n CT it(N c ), we rely on the basic observations and methods of Ref. |^1|. For large sizes of the 



box (AhL ^> 1 with Ah ~ 1 Gev), heavy excitations are exponentially suppressed in the 
partition function, which is then dominated by the lightest states, the pseudo-Goldstone 
bosons of SB^S. This leads to an effective description in terms of the Chiral Perturbation 



Theory (%PT) |2^j23|1 , and it can be matched with QCD, yielding the desired information 
concerning the infrared properties of the Dirac spectrum. Moreover, the effective Lagrangian 
is identical to its infinite- volume counterpart, provided that periodic boundary conditions 
are used 124] . 



If Nf lies far below n cr i t (N c ), the quark condensate is large and behaves at large (but 
finite) volumes according to the asymptotic behaviour derived by Leutwyler and Smilga 
PHI , using Standard |p2]| . Above n CT i t (N c ), the quark condensate vanishes, and the 

previous analysis cannot be applied. However, if chiral symmetry is still partially broken, 
the matching with remains possible and it leads to a clear-cut change in the large- 
volume behaviour of Ck : expressed through their inverse moments, the average behaviour 
of the lowest eigenvalues for L —>■ oo should turn from 1/L 4 ((qq) 7^ 0) into 1/L 2 ((qq) = 0) 
||13|| . When we approach the critical point with Nf near but under n crit (N c ), significant 
discrepancies from the asymptotic limit L —>■ 00 could be seen for large but finite boxes. 



The latter should then be analyzed using the framework of Generalized |23],|35|] . We have 
clearly in mind the possibility to use unquenched lattice simulations, varying Nf and (finite) 
lattice size L to eventually detect chiral phase transitions, through the volume dependence 
of inverse moments (|1.4| ). 

This article is organized as follows. In Sec. 2, we briefly review features of Euclidean 
QCD and of the effective theory on a torus. Sec. 3 explains how both theories are matched 
to derive the original form of Leutwyler- Smilga sum rules below n cr i t (N c ), before analyzing 
how they are modified in the phase where the quark condensate vanishes. In Sec. 4, we 
discuss the approach to the critical point, where a competition between a small quark con- 
densate and higher order contributions leads to sizeable computable finite-volume effects. 
Sec. 5 is devoted to the computation of the next-to-leading-order corrections to the sum 
rules. We discuss in Sec. 6 how to obtain from the inverse moments an unambiguous signal 
indicating that Nf approaches n cr i t (N c ), and we discuss the interest of lattice simulations in 
this framework. Sec. 7 summarizes the main results of this work. 
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II. SMALL MASS AND LARGE VOLUME EXPANSION OF THE PARTITION 

FUNCTION 



A. Euclidean QCD on a torus 

The Euclidean^ QCD Lagrangian for Nf light quarks reads: 

£(JV/) = T~^ G ^ ~ l9v ~ + Wq, (2-i) 



with the winding number density: 



1 



= ^G%(x)G%(x), (2.2) 
and the vacuum angle 9 The quark mass matrix M is of the form: 

M = 1(1 - 75 )M + 1(1 + 75 )Mt, (2.3) 

where M is a Nf x iVj complex matrix, diagonal in a suitable quark basis with positive real 
eigenvalues. 

We consider the partition function of this Euclidean theory in a finite box L x L x L x L, 
large enough to neglect safely the heavy quarks: 

Z\N } ) =C J [dG] j [di>][d^} exp (- j d 4 x C (N ^ , (2.4) 

where C is a normalization constant, which may depend on the volume, but not on the mass 
matrix. 

We impose boundary conditions on the fields, by viewing the box as a torus and identi- 
fying and Xp + n^L (with n M integers): the gluon fields have to be periodic and the quark 
fields antiperiodic in the four directions, up to a gauge transformation. The gauge fields 
are classified with respect to their winding number v = j v dx v{x), which is a topologically 
invariant integer (related to the gauge transformation defining the periodicity of the fields 
on the torus). The index theorem asserts that v is the difference between the number of 
left-handed and right-handed Dirac eigenvectors with a vanishing eigenvalue. 



The Dirac eigenvalues satisfy a uniform bound [11]: 



i/d 

\\n[G}\<C—=UJ n . (2.5) 

This bound means essentially that an external gauge field lowers the eigenvalues of the free 
field theory ||10|| . It involves a coefficient C, depending only on the geometry of the space- 



time manifold, but neither on G, n or V = L d . The partition function Z can be decomposed 
in Fourier modes over the winding number: 



2 In this paper, all the expressions are written in the Euclidean metric, unless explicitly stated. 
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v=—oo 

Each projection of positive winding number v is: 



(2.6) 



Z v (N f ) = C f [dG]e~ SYM[G] det(-ifl + M) 

J V 

= C f[dG]e- SYM[G] (det f MYY[ 



det/ (A 2 + MMt) 



(2.7) 
(2.8) 



/ denotes the integration over the set of the gluonic configurations with a fixed winding 
number u, and SVm is the pure gluonic action, det / M is the determinant of the Nf x Nf 
quark mass matrix (it is replaced by (det/ M^)~ u for v < 0). The primed product includes 
only the strictly positive eigenvalues: its denominator involves the Vafa-Witten bound u n 
of Eq. ( |2.5| ) and a reference mass scale \i larger than any light quark mass. It represents a 
convenient normalization of the determinant, such that each factor of the primed product 
is lower than 1 when the quark masses tend to zero. This normalization does not affect any 
observable. We check that the quark mass matrix and the vacuum angle arise in the partition 
function through the product M exp(iO/Nf), consistently with the anomalous Ward identity 
for the singlet axial- vector current. 

The partition function for a fixed positive winding number is: 



Z v (N f )=C f[dG]e~ SYM[G] (det MY 



(2.9) 



n A - 



1 . . „2 I eX P 



where () denotes the trace over flavours. Provided that the partition function is regularized, 
we can expand the logarithm for small masses (compared to the size of the box): 



Z v {N f ) =C f [dG\e 



«?">•( IT;st„, 



(2.10) 



x exp 



C' u (det M) v 



(M^M)a 2 - -((M+M) 2 )^ + 0(M 6 ) 
1 + (MtM) ((# - \{{M ] MY) ((a 4 »^' ; 



+ i(MtM) 2 «(a 2 ) 2 »^ ) + 0(M 6 ) 



(2.11) 



The inverse moments are defined for each gluonic configuration as: = 52 n 1/A£. The 
normalization factor C' u is independent of the quark mass matrix: 



C' U =C [dG]e 

J V 



— Sym[G] 



n 



(2.12) 
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The average over gluonic configurations with a given winding number is defined by 



where the denominator is a normalization factor, ((1))^ = 1. In Eq. (|2.11| ), this average is 
applied to inverse moments <Jk that are particularly sensitive to the infrared tail of the Dirac 
spectrum. On the other hand, (Q)^ f ^ includes a product over eigenvalues, which should 
suppress the statistical weight of the lowest eigenvalues when the number Nf of massless 
flavours increases. The averaged inverse moments in the exponential of Eq. (|2.11| ) could 
therefore exhibit a strong dependence on Nf. 

( |2.10|) contains several sources of divergences. Let us first consider the gluonic configu- 
ration as a fixed external field. In the fermion sector, sums over the Dirac spectrum may 
diverge because of its ultraviolet tail. For A — ► oo, the number of eigenvalues in [A, A + AA] 
is given by the free theory: 

An = -4f|A| 3 AA. (2.14) 

47T 2 ' ' V ' 

The expected ultraviolet divergences of the inverse moments have therefore to be subtracted. 
We can write: 

a 2 = a 2 + D ( 2 Nf \ a 4 = a 4 + Df f \ (2.15) 

where the divergent part is included in D, and a is finite. For instance, we can choose an 
ultraviolet cutoff A and define the integer K such that ojk = A. The regularized inverse 
moments then read: 



K 



n=l 



and the divergent parts behave (at the leading order of the volume) like: 

D ( 2 Nf) ~VA 2 , Df^-VlnA. (2.17) 

These short-distance contributions are the same for all winding-number sectors. If we per- 
form this splitting in Eq. ( 2.10| ), we obtain the regularized partition function Z v involving 



the inverse moments a, multiplied by an exponential factor with divergent counterterms 
which contribute only to the vacuum energy: 



Z v (Nf)=Z v (N f )exp 



D ( 2 f \M^M) - -D[ Nf) ((M^M) 2 ) 



(2.18) 



Secondly, the product over the eigenvalues in the fermion determinant of Eq. ( |2.13| ) 
needs a regularization already for a fixed gluonic configuration. Nevertheless, for observables 
dominated by the lowest Dirac eigenvalues, we expect less sensitivity to the ultraviolet tail 
of the determinant. If we split the product over eigenvalues into ultraviolet and infrared 
parts |||32|]: 
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K f \2 \ty 

A = A IR A UV , A m = TT ( -TT-2 ) > ( 2 - 19 ) 



n=l 



we can expect the gluonic average of the inverse moments to depend essentially on Am, with 
a weak sensitivity on A. 

Up to now, the gauge configuration was viewed as an external field, but the integration 
over the gluonic fields leads to a third series of divergences. Fortunately, their regularization 
is rather disconnected from the fermion sector [^U] (for instance, the cut-off may be chosen 
independently of A). For the purpose of this paper, it is sufficient to stick to a multiplicative 
renormalization of the mass matrix and the Dirac eigenvalues, 

M -> Z m M, A n -> Z m X n , (2.20) 

inducing a multiplicative renormalization for ((o*.))^ ■ We shall only consider homogeneous 
quantities, like ratios of inverse moments with the same degree of homogeneity in A: the 
problem of the renormalization in the gluonic sector is therefore discarded in the rest of this 
article. 



B. Effective Lagrangian 



For large volumes, the massive states are exponentially suppressed. The partition func- 
tion is therefore dominated by the Nj — 1 pseudo-Goldstone bosons resulting from the 
Spontaneous Breakdown of Chiral Symmetry and described at low energies by the Chiral 
Perturbation Theory (%PT). The effective Lagrangian for Goldstone bosons is written as a 
double expansion in powers of the momenta p and of the quark masses m: 



C 



eff 



5> 



(k,i), 



(2.21) 



kj 



where C(k,i) gathers all terms contributing like p k m l . In Euclidean QCD, it has been shown 
that, on a large torus, the low energy constants in £ e fr are not affected by finite-size effects 
21]. 

If U{x) G SXJ(Nf) collects the Goldstone fields, the partition function is: 



Z\N f ) = j [dU] exp - J d 4 x C^ f) (U, dU, Me i9/N f) . 



In this framework, the projection on a given winding number yields [[H 

dO 



1 



-ivQ 



[dU] exp 



— / [dU]{detU) v exp 

Z7l 



[ d 4 x C%<\u,Me»' N t] 
Jv 



- / d 4 xC$\u,M) 
v 



(2.22) 

(2.23) 
(2.24) 



with U(x) = U(x)exp(—id/Nf). The path integral over SXJ(Nf) for the partition function 
Z ends up with an integral over U(Nf) for Z v . Because of the invariance properties of the 
measures [dU] and [dU], we have for any Vi, V2 G XJ(Nf): 
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Z v {N f \V 1 MV 2 ) = (detVM)" Z v (N f \M). (2.25) 

The low-energy constants in C e g are volume-independent and iV^-dependent order param- 
eters. In particular, a partial restoration of chiral symmetry would make some of them 
vanish. Since the relative size of these order parameters vary with Nf, the organization of 
the double expansion ( |2.2ip depends on the phase in which the theory is considered. 

1. If the number of light flavours Nf is fixed below n cr i t (N c ), the quark condensate Y,(Nf) 
is the order parameter that dominates the description of SB^S for sufficiently small quark 
masses (or sufficiently large volumes). The leading order of the effective Lagrangian involves 
only a kinetic term and a term linear in the quark mass matrix: 

c (N f ) = ^(N^rf^u) _ ±X(N f )(rfM + MiU). (2.26) 

F is the decay constant of the Goldstone bosons and S(iV/) is the quark condensate, in- 
troduced in Sec. | in Eqs. ( |1.3|) and (1.2). The expansion of the effective Lagrangian is 



organized in this case through the standard power counting JZIJ: d ~ p, M ~ p 2 , so that 
the next-to- leading order is 0(p 4 ). 

2. On the other hand, for Nf > n cr i t (N c ), the quark condensate vanishes and we cannot 
rely on the previous description anymore. In this case, the leading-order Lagrangian is the 
sum of the kinetic term, C(2,a), an d of a term quadratic in the quark masses, £(0,2) : 

= ^p-(d»U%U), (2.27) 

C m = -7 [A(N f )((U^M) 2 + (M^U) 2 ) + Z s (N f )(U^M + M^U) 2 (2.28) 
+Zp{N f )(U ] M - M ] U) 2 + H{N f )(M*M)] . 

£(0,2) appears in the standard 0(p 4 ) Lagrangian at the next-to-leading order, and the low- 



energy constants Z$, Zp, A and Tt correspond respectively to L 6 , L 7 , L$ and H 2 of Ref. J22 
In this phase, the counting used to perform the expansion at higher orders is modified [33]: 
d ~ M ~ p. 

In the generic case Nf > 3 (the case of two flavours is commented in App. |B|), A, 
Zs and Zp are order parameters of SB%S. They are related to the low-energy behaviour 
of two-point correlators of the scalar and pseudoscalar densities^ S a (x) = ip(x)t a ip(x) and 
P a (x) = ^(x)i a i75^(x), where {t a } are flavour matrices. A stems from (S a Sb — P a Pb)- Z s 
is given by the correlator (SoSo^ — S a Sb), and Z P by (P -Po^a6 — P a Pb)'- Z s and Z P violate 
the Zweig rule in the scalar and pseudoscalar channels respectively. 

H, is a high-energy counterterm, which is not an order parameter and cannot be measured 
in low-energy processes. Other similar counterterms arise at higher orders: they involve 
only the quark mass matrix M, but not the Goldstone boson fields U. These counterterms 
are needed to subtract short-distance singularities in QCD correlation functions of quark 



3 Notice that contrary to the convention used in Refs. (2*| and |3l|, the decay constant F 2 is not 
factorized in £(0,2) : A, Zg and Zp carry the dimension (mass) 2 . 
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currents. Their general structure is dictated by the chiral symmetry, and it is reproduced 
by the high-energy counterterms on the level of the effective Lagrangian. 

3. For Nf just below the critical point n crit (N c ), we expect a small (but non- vanishing) 
condensate and a large Zweig-rule violation in the scalar sector |J. Linear and quadratic 
mass terms in the effective Lagrangian may be of comparable size. To take into account this 
possibility, we include both of them in the leading order of the Lagrangian: 

Cf f) = i [F 2 (N f )(d^d^U) - 2E{N f ){tfM + M ] U) (2.29) 

-A{N f )({U ] M) 2 + {M ] Uf) - Z s {N f )(U ] M + M ] U) 2 
-Zp(N f )(U^M - M^U) 2 - H(N/)(M*M)] . 

This Lagrangian can be actually viewed as the lowest order of another systematic expansion 



scheme, defined by the generalized chiral counting [23] : d ~ M ~ B ~ 0(p). In this case, 
the next-to-leading order counts as 0(p 3 ). 

The standard and generalized counting rules are only two different ways of expanding 
the same effective Lagrangian: 

£ eff = £ 2 + £ 4 + . . . = £2 + £3 + . . . (2.30) 

At a given order in p, Generalized includes terms relegated by Standard to 

higher orders. At the lowest order, Eq. fl2.29|) can be applied even if the quark condensate 
dominates. On the other hand, the Standard xPT becomes inaccurate in the vicinity of the 
critical point where S ~ 0, whereas Generalized %PT may be more appropriate to describe 
the transition. 



III. LEADING LARGE- VOLUME BEHAVIOUR OF THE INVERSE MOMENTS 

A. Matching QCD and the effective theory 

If we analyze perturbatively the partition function ( [2.22j ), the only difference from the 
case of an infinite volume lies in the meson propagator, because of the periodic boundary 
conditions: 

1 pip* 
p 

where p^ = lisn^j L, with n M integers. The contribution of the mode p = in this propagator 
blows up when pions become massless [^7[] . Graphs containing such zero modes will diverge 
in the chiral limit, whereas the non-zero modes are suppressed in the large- volume limit: 
the fluctuations of the zero modes are not Gaussian and cannot be treated perturbatively. 
To cope with them, we split the Goldstone boson fields in two unitary matrices: U(x) = 
U Ui(x), where the constant factor U describes the zero modes and U\(x) the remaining 
non-zero modes. 

In a first approximation, the Gaussian fluctuations of JJ\ can be neglected and the path 
integral in Z reduces to a group integral over constant SU(Aj) matrices: 
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Z(N f ) = V f [dU ] exp \-V£ l £ f \Uo,Mexp(iQ/Nf) 

JSV{N f ) L 



(3.2) 



where [dUo] is the Haar measure over the group, and T> a normalization constant, indepen- 
dent of the mass. The projection on a topological sector Q2.24|) becomes: 



z„m = -v 



[dU ]{detU ) v exp -VC { C ^\U 01 M) 



V(N f ) 



(3.3) 



To simplify the notations, we replace Uq by U in the calculations at the leading order of Z v . 
In addition, the iV/-dependence of the low-energy constants will not be explicitly denoted 
from now on, unless its presence is mandatory for understanding. 

We want to expand Z v with respect to the size of the box and to the quark mass matrix. 



Actually, Eq. (|3.2|) tells us how to organize this from the expansion of C e ^. At the leading 
order, the partition function will depend on a simple scaling variable X = ML K . Below 
n CT it(N c ), we have k = 4 (c.f. Eq. ( |2.26| )), whereas the phase with a vanishing condensate 
2 (c.f. Eq. (|2.29|) ). For small X, the expansion of Z u reads: 



'Cjl 

yields k - 



: AC(det Xy\l + a u (X^X) 

+K(X^X) 2 + c v {{X^Xf) + 0(X% 



(3.4) 



where the coefficients N v , a u , b u , c u do not depend on M. This expansion is valid for v > 0: 
for a negative u, (detX^)^ arises instead of (detX)^. The calculations are very similar in 
both cases, and our future results can be translated for any winding number by writing \u\ 
instead of v. 

The QCD partition function was expanded as a polynomial in the quark masses in 
Eq. ( pTTlD , leading to: 



C' u L~ KvN f {detfX) 
1 



(3.5) 



L 2 
1 



(^) 2 «(- 2 ) 2 »r /) +o(x 6 ) 



2L iK 



By identifying the same powers of X in both expansions, we obtain relations between param- 
eters of the effective Lagrangian and the leading large- volume behaviour of inverse moments. 

When we compare Eqs. (|3.4j ) and (|3~5|), we have to take into account the divergences of 
the inverse moments a^, as stressed in Eq. (|2.18| ): 



Z u exp 



Z* PT . 



(3.6) 



These counterterms, built from traces of the quark mass matrix, are also present in the 
%PT expression of the partition function. Therefore, the divergent behaviour of the inverse 
moments (e.g D 2 for a 2 ) is tracked by counterterms in the Lagrangian (in this case, T~C). 
Divergence-free sum rules are found by considering linear combinations where the related 
%PT counterterms cancel. 
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B. Nj <C n CT i t (N c ): Leutwyler-Smilga sum rules 



This case has already been treated with great details in Ref. [21]. We briefly review 
the main steps of the derivation of Leutwyler-Smilga sum rules for the reader's convenience. 
Eq. ( |3.3|) yields at the leading order: 

1 



2tt 



-V 



[dU}(det U) u exp 



V(N f ) 



{U ] M + M^U) 



(3.7) 



VMT, is the only parameter of the group integral, and the scaling variable is X = ML 4 
(k = 4). In the general case of an arbitrary matrix M, a formula for the integral (|3.7| ) is 
discussed in Ref. [ Z8f| . For our present purpose, it is however sufficient to follow the original 
method described in Ref. \21\ to expand Eq. (3.7) in powers of M. We obtain the expansion 
coefficients a u , b u , ... through two derivative operators, applied on both expressions of Z v : 



(3.8) 



the group integral ( |3.7| ) and the X-expansion (|3.4j) . The latter gives: 
d d 



E 



d x ad x* z ^^ detx y 



N f K 



i v + (X^X) [(N f K + l)b v + (Nf + K)c) + 0(X A ) | 



and 



^2(t a t b t c t d ) f v Q t a d v , z v = K(detxy 

abed 



dX n dX*dX c dX*' 



(3.9) 



X- 



NfK 



{(N f + K)b u + (NfK + l)c v + 0(X 2 )} 



where K = Nf + u, and X a are the coordinates of X on {t a } (a = . . . Nj — 1), which is a 
complete set of Hermitian matrices (see App. [A]). 

The same derivative operators are applied on the group integral ( ft.7|) : 



E 



abed 



_d d_ 

OX a dX*' 



d d d d 
dX a dX* b dX c dX* d 



— NfJ^Z: 



256 



(3.10) 
(3.11) 



Once Z v is replaced by its X-expansion ( p.4|) on the right hand-side of Eqs. ( |3.10|) and ( p.ll|) , 
these equations yield polynomials in X, which are identified with Eqs. ( |3.8| ) and ( |3.9j ) order 
by order in powers of X. We get thus a u , and a linear system of two equations for b u and 
c v . 

Once a„, b v and c v computed, the comparison of Eqs. ( |3.4| ) and ( |3 .5| ) leads to the 
Leutwyler-Smilga sum rules: 

, jVZ(Nf)} 2 



(N f ) 



\{°2?))l Nl) 16{K2 _ Vj 



2K 



AK 
[VX(Nf)} 4 



16K(K 2 



(3.12) 
(3.13) 
(3.14) 
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Because of K = Nf + \u\, the sum rules (|3.12|) -( |3.14|) depend explicitly on the number 
of flavours, but there is another (implicit and unknown) dependence stemming from the 
quark condensate E(iVj). No divergent counterterm is explicitly present: these sum rules 
are derived from the leading order Lagrangian in Standard xPT, and they show only an 
asymptotic behaviour, valid for V — > oo. For instance, 02 and (cr 2 ) 2 contain divergent 
subleading terms[]. 

C. Nf > n cr [ t (N c ): the phase with a vanishing quark condensate 

For Nf > n cr i t (N c ), the integral defining Z u in terms of the effective Lagrangian Q3.3| ) 
involves quadratic mass terms at the leading order: 



z AN f ) = ^-V [ [dU](detU) v (3.15) 

^ Jv(N f ) 

{A{{U ] Mf + (M^U) 2 ) + Z S (U ] M + M^U) 2 



'U(JV» 

'V 



x cxp 
• 4 

+Z P (IF~M - M ] U) 2 + H(M^M) } 

The scaling variable is now X = ML 2 (k = 2). The counterterm TC has the same structure 
as the divergent term D2 due to 02 in Eq. ( |3.6| ). To eliminate this divergence, it is natural 
to introduce the //-dependent fluctuation a 2 = 02 ~ {{ a 2)) u N ^ ■ The subtraction of this 
quadratic divergence leads to the loss of a single sum rule, for instance (((72))^^. For the 
other topological sectors, we can indeed write sum rules concerning ((c 2 ))^ V/ '' — ((o^o^) 
since the (short-distance) divergence due to TC is insensitive to the (topological) winding 
number. 

Because of chirality, the integral ( |3.15 ) vanishes unless the same power of U and TP arises. 



The determinant (det U) v counts as the vNf-ih power of U, whereas the exponential involves 
only the square of W. Therefore, the phase of an odd Nf > n crit (N c ) discriminates between 
the topological sectors: the odd-v sectors are suppressed in the large-volume limit compared 
to the even winding numbers (this discrimination does not occur for an even number of 
flavours). As a matter of fact, the symmetry M — > —M is equivalent to 9 — » 9 + nNf. From 



the Fourier decomposition (|2.6|) , we can directly check that the odd topological sectors have 
a vanishing partition function at the leading order, provided that Nf is odd. Of course, 
higher orders of the effective Lagrangian (for instance £3) contribute to the odd topological 
sectors, giving finally rise for Z v to a different volume dependence from the even winding 
numbers. 

In the topologically trivial sector v — 0, (det U) u disappears from the group integral and 
the exponential in (p.15 ) can be directly expanded in powers of X and integrated over U(N). 



4 For this reason, the formulae ( 3.12| ) and (3.13) should be applied to finite volumes with great 
care. 
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Using App. |A], the computation of the lowest powers in the X-expansion is straightforward, 
leading to the sum rules: 

«<^ 2 »° = 16iV^|-l) < 3 ' 16 ' 
x [M2N'f + 1)(Z| + Z%) - SZ S Z P - SN f A(Z s + Z P ) + 4N]A\ 

^ = W$f=T) (3 - 17) 

x [12(Z 2 S + Z 2 P ) - 8Z S Z P - 8N f A(Z s + Z P ) + AA 2 }. 

As emphasized in the previous section, these sum rules depend on the number of massless 
flavours in an explicit way, but also implicitly through the iVj-dependent order parameters 
A, Zs and Zp. 

These sum rules predict a different large-volume behaviour from the Leutwyler-Smilga 
sum rules ( |3.12| )-( |3~T4 ). This agrees with our general expectation concerning the large- 



volume dependence of the (suitably averaged) small Dirac eigenvalues [[y| . The eigenvalues 
accumulating like 1/L 4 contribute to SB^S and to the quark condensate. Correspondingly, 
for Nf < n CT i t (N c ), the asymptotic behaviour of the sum rules is: 



V\ >:il))~V\ (3.18) 





On the other hand, the l/L 2 -eigenvalues do not contribute to the quark condensate, but 
may still contribute to SB%S in the phase above n cr it(N c ), through a non- vanishing value of 
F 2 (Nf). Indeed, Eqs. ( |3.16|) and Q3.17| ) predict in this phase an infinite- volume limit of V 2 
for (((a- 2 ) 2 ))o and ((cr 4 )) , as expected. 



IV. THE APPROACH TO THE CRITICAL POINT 
A. Leading large-volume behaviour 

We want now to study the intermediate case, where the linear and the quadratic mass 
terms in the effective Lagrangian may compete for some range of volumes. To understand 
which results can be expected, it is instructive to consider first xPT in an infinite volume 
and to imagine that we let the quark masses vary. If the quark condensate is (even slightly) 
different from zero, we can always find sufficiently small quark masses for which the linear 
mass term is dominant. When the quarks become massive, the corrections due to the 
quadratic mass terms may become discernible and even preponderant, provided that the 
quark condensate is not too large to hide their effects. 

In this paper, we work in a box with a fixed large volume, and M 2 is counted as 0(1/L 4 ). 
The variation of the quark masses is therefore translated into a change of the volume. For 
Nf < n crit (N c ), the Leutwyler-Smilga sum rules derived in S%PT should correctly describe 
the volume-dependence of the inverse moments when L tends to infinity. However, close to 
the critical point and for a given value of the volume, the quark condensate need not be 
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large enough to make £ 2 , Eq. ( |2.26|) , dominate. This could lead to significant deviations 
from the asymptotic limit even for large volumes. 

Hence, the leading order of the Lagrangian is £2, Eq. ( |2.29| ), and Z v reads: 



Z v =^-v! [dU}(detU) u (4.1) 

^ J\J(N f ) 

x exp — {2S(f/ t M + M*U) + A((U^M) 2 + (M^U) 2 ) 

+Z S (U ] M + M ] U) 2 + Z P (WM - M ] U) 2 + H(M^M)} 

X = ML 2 remains the scaling parameter for the mass, and SL 2 is the expansion variable 
for the quark condensate. This organizes the expansion through the power counting E ~ 
M ~ 1/L 2 , similar to G^PT. We shall therefore consider the theory for volumes and masses 
so that X and £L 2 are of order 1. 



In order to evaluate ( |4.1|) , it is convenient to define the group integral I u for arbitrary 



complex numbers (6, b, z, z, y, a, a): 



I u = / [dU] (det U) v (4.2) 

JlJ(N f ) 

x exp[b(Xtf) + b(X^U) + z(XU^} 2 + z(X ] U) 2 

+y(XU^) (X?U) + a((XW) 2 ) + a((X^U) 2 )], 

The partition function at a fixed winding number reads: 

Z u = ^Vexp[h (X^X}}I u (b ,b ,z ,z } y ,a } a°;X), (4.3) 
where I u is calculated with the values: 

b° = b° = ^L 2 S, (4.4) 

z° = z° = \(Z S + Z P ), y° = \{Z S -Z P ), (4.5) 

a = d° = \A, h° = \n. (4.6) 
4 4 

l v is a polynomial in (b, b, z, z, y, a, a), and its derivatives are not independent: 



d 2 i v di v d 2 i u di v d 2 h di v 



db 2 dz db 2 dz dbdb dy 



(4.7) 



We expand this integral in powers of X, with coefficients that are independent of the quark 
mass matrix: 

I v = (det X) v [a v + (3 V (X^X) + ^(X^X) 2 + 5 U ({X^X) 2 ) (4.8) 
+e u {X^Xf + ^((XtX) 2 )(XtX) + K u ((X^Xf) + 0(X% 
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We identify the same powers of X in the expression of Z v in terms of averaged inverse 
moments (|3.5| ) and in its expression at the leading order of the effective Lagrangian, obtained 
from Eqs. (|4.3|) and ( fO ). This leads to the sum rules: 



«M 2 » 



V 

2V 2 



0u 



+ h 



lv , , A, .h 2 

h n 

2 



a 



ft,. 



-2V 2 ^ 



\ I J— - I ^ 



ft,. 



ft,. 



V I 6^-6^^ + 21 — 



ft,. 



({W^-V* (-2^ + 2^ 



W 6 K h 



(4.9) 
(4.10) 
(4.11) 

(4.12) 

(4.13) 

(4.14) 
(4.15) 



If we know a u , /?„,. . . in terms of the low-energy constants of C 2 , Eqs. (|4.9|) - ([4.15 ) lead to the 
desired sum rules. The high-energy counterterm h, which reflects the ultraviolet divergence 
in 02, has to be eliminated. This can be obtained if we consider the fluctuation of a-i over a 
topological sector: <7 2 = o"2 — (( cr 2))[ / 7V/ \ as defined in Sec. |111 C. 



For the topologically trivial sector v = 0, the computation is very simple, following the 
same line as for the phase Nf > n crit (N c ). This leads to the expansion coefficients (for b = b, 
z = z, a = a): 



fto 
7o 



0o 
1 



1 



(y + b 2 



N f (N 2 - 1) 

x { N f[~2 + 2b2y + 2b2z + y 2 + 2zl + 2fl2 ] - 2a ^ 2 + 22 l} 



(4.16) 
(4.17) 



So 



N f (N 2 - 1) 

x + 26 2 ^ + 2b 2 y + y 2 + 2z 2 + 2a 2 ] + N f ■ 2a[b 2 + 2z] J 



(4.18) 



eo 



Nf(N 2 - l)(iV 2 -4) 
6 6 



(4.19) 



x <^ 6(iV; 



36 V3 4. 
+2(iV 2 + 2) [b 2 +y] a 2 -12N f 



— + o I — + 7 I +o [ z + -— + yz + I ^- + 
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«0 



N f (Nj - l)(JVj-4) 



x 



6 6 



-18N f 

36 V3 I 

-18iV/ [b 2 + y] a 2 + 12(N 2 f + 2) 



J + b2 (l + z ) +yz 



Nf(Nj - l)(iV|-4) 



24 



36 



£ ,V 
3 4 



2 i y i . 

2 + y + 3/2 J + 



2/ . 2 
6 



+A(N 2 f + 2) [6 2 + y] a 2 - 2AN f 



~6 



+ b 2 ' y 



+ z 



yz 



(4.20) 



(4.21) 



(4.22) 



Before focusing on the resulting sum rules for the topologically trivial sector v = 0, we 
sketch the general derivation of the expansion coefficients for an arbitrary winding number. 



B. Topologically non-trivial sectors: v ^ 

Let us begin with the leading coefficient a v . Independent of X, it can be computed for 
X = x ■ 1, where x is a complex number. a v is then given by the leading order of I v in x 
(without any power of x*), and it depends only on (6, z, a). As a matter of fact, a u (b, z, a) 
can be deduced from a v (b, 0, a) because of the relations between the derivatives ([4.71). The 
problem reduces to obtaining the leading order in x of the group integral: 

I" = I v {b, a; x ■ 1) = / [dU}(detUyexp[bx(U j )+ax 2 (U^ 2 )}. (4.23) 

JV{Nf) 



The Appendix |C 1| describes how a v (b, 0, a) is extracted from this integral, leading to the 
polynomial: 

a u (b,z = 0,a)= b vN !~ 2m a m c m , (4.24) 

m=0...uN f /2 

where {c m } are purely combinatorial coefficients. Using d 2 a u /db 2 = da u /dz, we obtain the 
general expression of a v : 

a u (b,z,a)= feW> (/ + y )! c m . (4.25) 

In the limit case of a vanishing quark condensate (6 = 0), we check that a v (and therefore I u ) 
vanishes if uNf is odd, in agreement with the parity discrimination discussed in Sec. [Ill C| . 

We obtain the next coefficients by applying the derivative operators of Eqs. ( |3.8| ) and 
( p.9[) on both representations of l v : the group integral (|4.2| ) and the X-expansion (|4.8|) . We 
already know the result of the latter from the phase Nf <C n CTit (N c ), studied in Sec [II [B[ 



18 



a " 

= (detX) 1 



(4.26) 



N f K 



P v + (XtX) [{N f K + 1) 7 , + (N } + K)5 U ] + 0(X 4 )} 



E, . d d d d 

{t a t b t c t d ) "7TTT - „„. TTTT" „„. i,y 



abed 



dx a dx; ox c dx* 



(4.27) 



(detX)"^- {(JV) + X) 7 , + (iV/K + + 0(X 2 )} . 
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The two-derivative operator, applied on the group integral ( |4.2| ) that defines I u , leads to: 

Y——i 



-^{y + 66) + 2aa(X t X) 

- N f d N f - d 

+ (N f zb + ab + -jby)— + {N f zb + ab + -y&l/)^ 



+ (2iV)zz + — } -y 2 + 2az + 2az) — 
2 oy 



L. 



(4.28) 



We can now replace I u by its X-expansion (|4.8| ), and identify the resulting polynomial in X 
with the right hand-side of Eq. (|4.26|) . When we identify the coefficients of X°, we obtain 
(3 U in terms of a v and its derivatives: 



a, 



da v 
~db' 



da 



a,. 



d 2 a u 



dz 



(4.29) 



The coefficients of (X^X) lead to an equality between a linear combination of 7^ and S v , 
and some derivatives of a u and (3 V (these derivatives can actually be rewritten only in terms 
of derivatives of a u , since we know how (3 V is related to a v ). 

We follow the same line with the four-derivative operator. Actually, when we apply the 
operator to the group integral (|4.2), we only need the lowest power of X, to compare it with 
Q4.27 ). Factors of higher degrees, similar to aa(X'X) in Eq. ( |4.28|) can be ignored, and we 
obtain: 



E, . d d d d 

{tatbtctd) --- ~ -- - — 1, 

abed 



dX a dX*dX c dX* d 



(4.30) 



-6 2 6 2 + b 2 (N f a + z) + b 2 (N f a + z) + 2bby 



+y 2 + 2zz + 2aa + 2N f (az + az) 



[bb(N f by + 2N f bz + 26a) + b(QN f yz + Aay + 2N 2 ay) 

d 

+b(2N f y 2 + AN 2 az + AN f zz + AN f aa + Aaz)] — 
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+- 



-b 2 y 2 + b 2 (2N f z 2 + Aaz) + bb(AN f yz + Aay) 
+y 2 (5N f z + (Nj + A)a) + ANfaz 2 

d 2 



+8N f aaz + AN f z 2 z + Aa 2 a 



J db 2 

I - d 3 1 d 4 

+ - [2N f byz 2 + N f by 2 z + aby 2 + Aabyz] — + - [N f z + 2a]zy 2 



4 



d 2 



+ ^ + 2 z] a 2 l + la%y dadb 



db 3 

1 2 2 d 3 

+ A ay dadb 2 



«% 4 



+ 0(X< 



We replace I v by its X-expansion on the right-hand side of this equation. We keep only 
the coefficient for X° and we compare it with Eq. (|4.27 ), to end up with a second equality 
relating a linear combination of r y u and 5 V to the derivatives of ot v . The resulting expressions 
are listed in App. |C 2\ but it seems difficult to handle them in general. 



C. Topologically trivial sector: v = 

From the expansion coefficients cto, /3q. . . of Sec. |IV A| , we get the sum rules for the inverse 
moments of degree 4 and 6. If we denote ( = VT? /A, S = Z$/A and P = Zp/A, the sum 
rules read: 

«<*>'»™ = Tum^T) < 4 - 3i » 

x {C 2 + ([A(2N 2 f + 1)5 - 4P - 4JV>] 
+ [4(2iV 2 + 1)(5 2 + P 2 ) - 85P - 8N f (S + P) + AN 2 ]} 

x{C 2 + C[125-4P-4X / ] 
+ [125 2 + 12P 2 - 8SP + 4 - 8N f S - 8N f P}} 

II WW - V3A3 (A QQ"\ 

(K))o " 32^(^-1X^-4) (433) 
x {C 3 + C 2 [305-6P-6iV / ] 
+C[1805 2 + 36P 2 - 12SP - 72N f S + 6(iVf + 2)] 
+ [1205 3 - 120P 3 + 725P 2 - 725 2 P - 72N f S 2 
+72N f P 2 + 12(iV 2 + 2)5 - 12(iV 2 + 2)P]} 

*«**™ = 16^(^1(^-4) (434) 

x {C 3 + C[2(2N 2 f + 7)5 - 6P - 6iV/] 
+C[36(iVj + 1)5 2 + A(Nj + 5)P 2 - 8(N 2 f + 5)5P 

-8N f (N 2 f + 5)5 + 4(2JVj + 1)] 
+ [24(iVj + 1)5 3 - 2A(N 2 + 1)P 3 + 8(N 2 + 5)5P 2 
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-8(N 2 f + 5)5 2 P - 8N f (N 2 f + 5)5 2 + 8N f (N 2 f + 5)P 2 
+8{2N 2 f + 1)5 - 8(2N 2 f + 1)P]} 

«<^ 3 »™ = mfiNf%^) < 435 > 

x {C 3 + C[Q(N 2 f + 1)5 - 6P - QNf] 
+C[6(2iV) -N 2 f + 2)S 2 + 6(iVj + 2)P 2 

-12(JVj + 2)5P - l2N f (Nj + 2)5 + 9iVj] 
+ [A(2Nf -Nj + 2)5 3 - A(2Nj - N 2 f + 2)P 3 

+ 12{N 2 + 2)5P 2 - 12{N 2 + 2)5 2 P 

-12N f (N] + 2)5 2 + 12N f (N 2 f + 2)P 2 

+ 18N 2 f S -18N 2 P}} 

The dependence on the number of massless flavours is not limited to the polynomials in Nf 
explicitly present in the previous formulae, since S, Zg, Zp and A are unknown functions 
of Nf (this dependence is here omitted for typographical convenience). The singularities for 
Nf = 1 (for 1 / A 4 - and 1 / A 6 -moments) and Nf = 2 (for 1 / A 6 -moments) arise because some of 
the coefficients a u , (3 V . . .in ( |4.8| ) are not independent in these cases, and we can only write 
(singularity-free) sum rules for differences between inverse moments of the same degree, e.g. 
(a 2 ) 2 - a 4 for N f = 1. 

Notice that the scaling volume parameter ( = VY? / A and the ratios 5 = Z$/A and 
P = Zp/ A are invariant under the QCD renormalization group. This invariance occurs also 
for ratios of inverse moments with the same degree of homogeneity in A: 

„ ((<74))o q (((^2) 3 ))o rp ((^ 2 ))p ((^ 4 ))q /2 u „ v 

(((^ 2 ) 2 ))o' K»o ' (Mo ' «a 6 » ' 1 J 

We can plot (Figs, [l] and §) the variation of R as a function of the volume, measured in 
physical units F~ A (P 7r =92.4 MeV). The scaling parameter is ( = (F^V)/ (16Lg), where the 
dimensionless parameter L 8 denotes (P 4 ^4.)/(16S 2 ) (for Nf = 3, it essentially corresponds 
to the S^PT low-energy constant L 8 of Ref. ||22||)p|. A variation of the condensate means a 
variation of Lg, and consequently a redefinition of the units used to measure the volume: this 
reduces to a simple shift of the curve (to the right if £ decreases, to the left if it increases). 

The infinite- volume limit reproduces the Leutwyler-Smilga sum rules (R —>■ Nf). On 
the other hand, since the scaling volume parameter is ( = VT? / A, the limit L —>■ corre- 
sponds mathematically to a vanishing condensate for the sum rules: we recover the results 
of Sec. |III Q . The sum rules Q4.31|) - (|4.35|) interpolate between these regimes. 



The ratios P, 5, T and U are not very sensitive to P (Zweig-rule violation in the pseu- 
doscalar channel) until we reach small volumes where large corrections stemming from higher 



5 In SxPT, the constant Lg depends on the renormalization scale fi. At fi = M p , it is estimated 
as Lg(Mp) = (0.9 ±0.3) • 10~ 3 (see for instance Ref. [^|). Close to the phase transition, L 8 should 
increase and become scale independent. 
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orders are expected. In the case of Nf = 3 flavours, the value P = —1/2 is privileged, be- 
cause it guarantees the validity of the Gell-Mann-Okubo formula, independently of the size 
of E. On the other hand, it may be interesting to notice that some ratios are affected by 
variations of S even at intermediate volumes. For instance, the dependence of the ratio S 
on S is plotted on Fig. |3] (we choose L 8 = 0.1, but other values of L$ can be obtained by a 
simple shift of the curve) . 

To simplify the analysis, it may be interesting to focus on linear combinations of the 
inverse moments in which the leading power of V cancels. These combinations therefore 
vanish in the limiting case of the Leutwyler-Smilga sum rules ( |3.12| )-( p.l4j ): 



Nf «(^ 2 ) 2 »r /) - Mo^ = ^& z s#V + [2Z 2 + 2Z 2 P + A}} (4.37) 



N f «^ 2 » ^ - 2 ((<T 6 ))f f) = mA V N , - 1} (4.38) 

x{2Z 5 SV 2 + [18Z| + 2Z 2 P - AZ s Zp - ANfZ s A + A 2 ]Y?V 
+[12(Z| - Z%) + AZ S Z P {Z S - Z P ) 

-AN f {Z 2 s - Z 2 P )A + 2{Z S - Z P )A 2 ]} 

Nf mX f) - 2 M ( " f) = sn/n]-!) (4 - 39) 

x{2Z 5 SV 2 + [Q(2N 2 f + l)Z 2 s + 2Z 2 P - AZ S Z P - AN f Z s A + A 2 ]Y?V 
+ [A(2N 2 + 1)(Z| - Z%) + AZ s Zp(Z s - Z P ) 

-AN f (Z 2 s - Z 2 P )A + 2(Z S - Z P )A 2 ]} 
Nf «(- 2 ) 3 »r /} - 3N f ((a 4 a 2 ))^ + ((a 6 ))^ (4.40) 
= ^-[3ZjVi: 2 + 2(Zl-Z P )} 

The large- volume behaviour of these combinations is particularly sensitive to the condensate 
S and to its fluctuation described by Z$ (Zweig-rule violation in the scalar channel). Both 
these parameters are precisely expected to be strongly affected by the vicinity of the critical 
point. 



D. Positivity conditions 

((T2) 2 , CT4 and <T6 are by definition positive, and their average over any topological sector 
should be positive as well. For Nf n crit (N c ), this positivity is trivially satisfied by the 
asymptotic behaviours predicted by the Leutwyler-Smilga sum rules (|3.12j) -( 3"7L4]) . 



When Nf is near (and below) n cvit (N c ), the volume dependence of the positive inverse 
moments is expressed through the sum rules of the previous section. They were derived 
at the leading order, for v — 0, and are functions of £, S — Z s /A and P = Z P /A. The 

positivity of (((<5"2) 2 ))o Vf \ ((^^o^ anc ^ ((°6))o P u ^ s therefore constraints on the low-energy 
constants of £ 2 - 

In the plane (S, P), it is instructive to draw the domain where each of these sum rules is 
positive for any value of ( = VY? j A: we demand the positivity of a polynomial of second or 
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third degree in (, whose coefficients are functions of S and P (and Nf). For a given number 
of flavours, this procedure excludes some values of (S, P), which constitute the hatched areas 
on Figs. ^-0. If (((o"2) 2 ))o iV/ ' > does not constrain S and P very much (Fig. |), the positivity 
of ((o"4))c/ Vf ' ) (Fig. H) and ((cq))^^ (Fig. leads to stronger conditions. If Nf increases, the 
excluded domains broaden, as shown on Fig. [7], compared to Fig. |5| If we suppose that 
Nf = 3 is already near the critical point n cr i t (N c ), and if we fix P = —1/2 from the Gell- 
Mann-Okubo formula, the positivity of ((o" 4 ))q 7V/ ^ yields the condition S > 1/6, explaining 
the zero in the plot of R on Fig. [l|, where the parameters S and P have been chosen on 
borderline of the positivity domain of ((o~4))^ f \ 

Obviously, these areas are obtained through the leading-order approximation to the sum 
rules: the border of these domains is altered by subleading corrections, which should become 
large for small volumes. Furthermore, the pseudo-Goldstone bosons do not dominate the 
partition function if the box becomes smaller than 1/Aqcd- To sum up, when we want 
the leading order of the sum rules to be positive for any (, we demand too much, and the 
resulting area is only an approximation of the really allowed domain in the plane (S, P). 

Furthermore, these positivity plots are only relevant for a number of flavour Nf ~ 

n CTit (N c ). Above the phase transition, (((o'2) 2 )) ( Q if \ ((c^o^ an d ((^e))^^ are still positive, 
but their large-volume behaviour is related in a different way to the low-energy constants 



of the effective Lagrangian, as described in Sec. [Ill C| . The positivity conditions stemming 
from the asymptotic behaviour of (((c*^) 2 ))^^ an d {(^e))^^ are satisfied for any S and P. 
The only non-trivial relation is due to the sum rule (|3.17|) for ((0-4))^^ and reads: 



(S + P — Nf) + 2(S — P) > Nj — 1. (4.41) 

To obtain this relation, we demand the infinite-volume limit of ((04)) to be positive. This 
limit is predicted by the sum rule ( 3.1 7]) at the leading order. The subleading corrections 



to this sum rule vanish as L — ► 00 and they do not affect ( 4.41|) . On the contrary to the 



previous positivity conditions obtained near the critical point, ( [4.41 ) is therefore exact for 
Nf > n clit (N c ). 



V. SUBLEADING CORRECTIONS 



This section is devoted to the next-to-leading contributions to the sum rules. In both 
phases, they behave as 1/L 2 compared to the leading order considered so far. 



A. N f < n crit (iV c ) 

The Leutwyler-Smilga sum rules were obtained at the leading order of the effective La- 
grangian £ 2 in the S%PT counting, restricted to the zero modes. The subleading corrections 
stem a priori from two sides: the non-zero modes (present already in £2), and the zero 
modes (beginning at the next-to- leading order £4). The first subleading corrections turn 
out to be of order 1/L 2 , and they come from the non-zero modes contributions to €2- They 
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can be expressed as a (volume-dependent) renormalization of the quark condensate^ in the 
sum rules (CT)- (CT) . 

The second type of subleading corrections arises from the zero-mode contribution to £4, 
quadratic in the quark mass matrix. This Lagrangian involves, among other terms, the 
counterterm (M^M) corresponding to the quadratic divergence of o 2 - Since the counting 
rule in this phase is M ~ 1/V, these quadratic terms are suppressed by a factor 1/L 4 in 
comparison with the linear term in £2. Consequently, they appear as next-to-next-to- leading 
order contributions and will not be discussed here. 

The non-zero modes arise in the decomposition of the Goldstone boson fields in Sec. [Ill A|: 



U(x) = UqU^x) = C/ exp li C{x)t a /F 



a=l 



2tt 

)~ exp \ l ~^~ n ' x 



(5.1) 



where is a four- vector whose components are integers (n ^ means ^ | rz^ | 7^ 0). The 
unitarity of U\(x) leads to = (0^)*. The fluctuations of the non-zero modes are small, 
leading to the counting rule £~0~<9~1/L. The leading contribution for the non-zero 
modes is <9 M £<9 M £ and comes from the kinetic term of £2. It is counted with the same power 
as the leading term of the zero modes (p.7|), but it can be directly integrated and becomes 
a simple contribution to the vacuum energy ||27| . At the leading order, the zero modes are 
actually the only relevant degrees of freedom. 

At the next-to-leading order, the corrections from the non-zero modes are due to the 
terms <9 2 £ 4 and M£ 2 . They are only suppressed by a factor 1/L 2 in comparison with the 
leading order L 4 M. The partition function ( |2.24|) up to the next-to-leading order is finally: 



Z v = V / [dOo](det u Q y e -^jj{u ,x) / FT d(f) a dm 

JV(Nf) 1 



(5.2) 



n>0,a 



x <>x i } 1 ~ £ M 

nj^0,a,b 



o \ 2 2 

2n \ n z 



T 4lP + 0{L- 



where the condition n > means: n > 0, or (no — 0,ni > 0), or (n = n\ = 0,n 2 > 0), or 
(no = ni = n 2 = 0, n 3 > 0). The (N 2 — 1) x (Nj — 1) matrix arises: 



Q 



ab 



2F 2 



((t a t b + t b t a )(U j M + M^U )). 



stands for the quartic term: 



2L 2 - 2 



7T 



■40 



3L 2 



E r n 4> b p 4> c q 4> d r n-(q-p) (t a t b t c t d ). 



(5.3) 



(5.4) 



abcd,npqr^O 



In the integral over the non-zero modes in Eq. (|5.2|) , these terms are suppressed by 1/L 2 
compared to the kinetic term. 



6 This result can be compared to the analysis performed in Ref. concerning the finite size-effects 
arising in the effective description of a spontaneously broken 0(A r )-symmetry. 
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We begin with the term T^, which involves neither the quark mass matrix nor the zero- 
mode matrix Uq. We can treat it perturbatively to perform an expansion in powers of 1/L, 
leading to: 



n d ^ d ^ 



(5.5) 



n>0,a 



x exp 



E< 



n \ 2 9 

T 



Sab + Q 



ab 




T44 + 



We should now apply Wick's theorem and contract the fields in T^. We would use the 
propagator stemming from the kinetic term and the "mass term" Q a b, where the latter is 
suppressed by 1/L 2 compared to the first. But we want only the first subleading correction 
due to the tadpoles arising in T^. Since this correction is already l/L 2 -suppressed compared 
to the leading order of the partition function, it can be calculated with propagators restricted 
to their momentum part (Q a b would induce l/L 4 -corrections). At the next-to- leading order, 
the contribution of T 4 ^ involves neither M, nor U (which are only present in Q a b)'- it is 
a global L-dependent term which can be factorized and eliminated by a redefinition of the 
normalization constant V . 

Hence, the 1 / L 2 -corrections are only due to the "mass term" Q ao of the non-zero modes. 
The partition function restricted to a given topological sector becomes: 



Z v 



V" I ldU }(detU yex P (-£ eff (Uo,X)) 

IU(N f ) 

I / \ - I > 1 2~ \ ~ 

x exp 



n 

n>0,a 



4 E 



n>0,a,b 



1 S. 



-n 
2 



ab 



Q 



ab 



0(L-^ 



The Gaussian integral over {</> n } can now be performed: 



(5.6) 



M JJexp 



;TrQ 



Af exp 



L 2 



4vr 2 



^2 n 2 



Tr Q 



where Af is an M-independent normalization factor. The trace over a, b = 1 
to: 



Tr Q 



N 2 — IE 



(MUi + UoM*). 



2N f F 2 

The integration over the non-zero modes ends up with the renormalization: 



X(N f )^E(N f )ll + g 



2N t 



9 



2n 2 F 2 L 2 ' 



E 



/ 1 



n- 



(5.7) 



N 2 - 1 leads 



(5.8) 



(5.9) 



If we include the first subleading corrections, the sum rules ( |3.12 )-( 5.14 ) remain therefore 
correct, provided that the parameters of the effective Lagrangian are renormalized, introduc- 
ing an additional l/L 2 -dependence related to the regularization scheme. In the dimensional 
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regularization introduced by Hasenfratz and Leutwyler fl3~0"f , the divergent sum v becomes 
— 47T 2 /?!, where f3\ is a "shape coefficient", related to the dimension and the geometry of the 
space-time. For a four-dimensional torus, (3\ = 0.1405 (see App. |B| for further comments). 

In this case, the first subleading corrections to Eqs. (|3.12| )- (|3.14| ) are summed up by the 
renor malizat ion : 

/ N 2 f -l R, \ 
E(JV» - E e (JV» = E(JV» (1 - -Lj- ■ ^- 2 j . (5.10) 

For instance, the relative correction (E — E c )/E remains smaller than a if the box size is 
greater than: 




L min = -J^—^ } (5.11) 



so that, for Nf = 3 flavours, the renormalization of E in the sum rules leads to a correction 
smaller than ten percent for box sizes larger than 1.9/ F (in the case of the dimensional 
regularization). 



B. Near the critical point 

As before, two sources of subleading corrections should be taken into account: the non- 
zero modes from £ 2 , Eq. ( 2.26|) , and the zero modes from the next-to-leading Lagrangian 



£ 3 @: 



Cf^ = I {X^Nf^d^d^fJiM^U + U ] M)) (5.12) 

+X(N f )(d^d^U)(M^U + U ] M) 
-7^i(A r / )((M t C/) 3 + {U ] Mf) 
-n 2 {N f )({M j U + U^M)M^M) 
-7l 3 (N f )(M*U - U ] M)((M ] U) 2 - {U ] M) 2 ) 
-TZ 4 (N f )((M^U) 2 + {U ] M) 2 )(M ] U + f/ f M) 
-n 5 (N f )(M^M)(M^U + U ] M) 
-TZ 6 (N f )(M^U - U ] M) 2 (M ] U + U^M) 
-TZ 7 (N f ){M ] U + £/ t M) 3 } . 

Since the counting rule is ML 2 ~ 1, both types of corrections are expected to con- 
tribute at the next-to-leading order 0(1/L 2 ), and could affect the previous quadratic or 
cubic volume-dependence of the sum rules. 

The non-zero modes are explicitly defined by ( |5.1|) . Like in the standard counting, their 
leading term in the effective Lagrangian is the kinetic term d^d^C,, which is counted as 
0(1/L A ). Its contribution (at the leading order) reduces to an overall constant, redefining 
the normalization of the partition function. 

The next-to-leading contributions from the non-zero modes are of the form B a M b d c ^ d , 
with 2a + 2b + c + d = 6, c and d even, and c < d. The possible terms are BMC, 2 , M 2 £ 2 and 
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<9 2 £ 4 from £2, and M<9 2 £ 2 from £3. At the next-to-leading order order, the path integral 
becomes: 



z v = v 



V(N f ) 



[dU }(detU o y eM-£ eff (U ,X)) 



II W d ^nT 
n>0,a 



(5.13) 



L 



x ex p - -v E w 



n>0,a,b 



T^ + O (- 



with the (iV| — 1) x (Nj — 1) matrices: 



Pab 
Qab 



2F 2 
1 

F 2 



;((t a t b + t b t a )(tf M + M^U )) 



(5.14) 



+A(t a UlMt b UlM + t a M^U t b M^U ) 
+Z s (t a (UtM - MiU ))(t b (U*M - M^U )) 

+Z P (t a (UtM + M^U ))(t b (uiM + M*U )) 



(5.15) 



The quartic term remains identical to its expression in the standard case (|5.4|) and it 
stems from the kinetic term of £2, whereas P is due to £3 and Q to the non-derivative 
part of £2- In Eq. (|5.13|) , the contributions of these three terms are suppressed by 1/L 2 , 
compared to the kinetic term: 7r 2 L 2 J2 n >o a ra2 |0n| 2 - 

For the same reasons as in the previous section, the integration of leads at this 
order to a term independent of M and Uq, which merely redefines the overall normalization 
constant T>. At the next-to-leading order, the partition function for a given winding number 
reads: 



v" I [du ){detu yeM-£ef f (Uo,x)) / TT drjm 

J n>0,a 



(5.16) 



x exp 



L 



n>0,a,b 

which yields after the integration over <fi: 



'n ( 1 



Af JJexp 



-Tr P- 



L 2 



4ir 2 n 



(5.17) 



M exp 



v n^0 



where M is an M- and [To-independent normalization factor. The traces are taken over the 




indices a,b — 1 . . . N 2 — 1: 
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Tr P 
Tr Q 



1 



F 2 L 2 2N 



2F 2 L i 



£L 2 



-{X + N f X)(XU 1 +U X^ 



(5.18) 



Nf-1 



+ (z s + 

2(Z S — Zj>) ^ ; : , .;- 



.4 

iV 



iV> 



(5.19) 



X 



/ 



(XC/ t )(Xtf7 ) - 2(Z S - Zp)(X^X) 



The integration over the non-zero modes ends up with a term of the same structure as £2, 
i.e. it renormalizes the parameters of the Lagrangian in the sum rules: 



A(N f ) 

Z s (N f ) 

Z P (N f ) 
H(N f ) 



E(N f )+g 



1 



-E(JV» + 



2m X 2 -l 



2JV) 

Z 5 (X / )+Z P (X / )- 

' r 4(Xv) _ ZpiNfl 
2 N f 
A(N f ) Zs{N f ) 
2 N f 
n(Nf)-g-2[Z s (N f )-Z P (Nf)} 



A(N f ) + g 

Z s {N f )+g 

Z P (N f )+g 



F 2 V 2N f 
A(N f ) 



[X(N f )+N f X(N f )}, 



N 



f 



with the sums to be regularized: 



2n 2 F 2 L 2 



n 2' 



u 



(5.20) 

(5.21) 

(5.22) 

(5.23) 
(5.24) 

(5.25) 



If we consider the dimensional regularization, we get g = —2f3i/F 2 L 2 and u = (see App. 0). 

With the counting £L 2 ~ ML 2 ~ 1, the first sub leading corrections stem also from the 
zero modes in £3: they contain therefore the low-energy constants TZi. If we consider the 
topological sector v — 0, the resulting corrections are quite simple to compute. When we 
expand Z (restricted to the zero modes) as a polynomial X, the integrals with different 
powers of U and U' vanish. In particular, the terms from £ 3 involve odd powers of the 
meson matrix and have to be combined with E(X'C7 + U'X). The resulting corrections are 
therefore TfIZi and are counted as 0(1/L 2 ). 

For v = the final form of the sum rules, including the first subleading corrections, is: 



V 2 A 2 



«(^2) 2 »0 



16X 2 (X 2 
V 2 A 2 



16N f (Nf - 1) K + s ? + a 4j> 



(5.26) 
(5.27) 



where s , is the leading term, already calculated in Sec. [LV C| , sj^ collects the terms from 
the zero modes in £ 3 , and s£ is due to the renormalization of £2 induced by the non-zero 
modes. The result is: 
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and 



s° 2 = ( 2 + ([4(2N 2 f + l)S-4P-4N f } (5.28) 
+[4(2JVj + 1)(5 2 + P 2 ) - 8SP - 8N f (S + P) + 4iVf], 

= "^[16^(^3 - K A ) + 8N f (N 2 f - 1)^5 - WN 2 f TZ 6 + 48iVfft 7 ], (5.29) 

s r 2 = JL {2(N 2 f - 1)C 2 + ([8N](N] - 1)5 - lQN 2 f P + 8N f ] 

+[-8(Nj - 1)(5 2 + P 2 ) - 16(3iV| + 1)5P 

+16iV>(iV^ + 1) (5 + P) -16N]]} (5.30) 

+ u ■ 4 ^ - 1 E( ^ + J^ } {< + 2 [(2JV? + 1)3 - P - N t ] } , 



s ° = ( 2 + ([125 - 4P - 4iV>] (5.31) 
+ [125 2 + 12P 2 - 85P - 8N f S - 8N f P + 4], 

sf = _[_8(JVJ - 1)^ 2 + 16JV>7Zg - - 16^ 6 + 48ft 7 ], (5.32) 

s l = jf { 2 ( N f ~ !)C 2 + C[8(N 2 - 1)5 - 8(iVj + 1)P - ANf(N 2 - 3)] 

+ [-8{N 2 - 1)(5 2 + P 2 ) - 16(JVj + 3)5P 

+32iV7(5 + P) -8(iVf + l)]} (5.33) 



It is worth commenting the above results: in the vicinity of the critical point, char- 
acterized by the counting EL 2 ~ ML 2 ~ 1, all the terms of the leading contribution s\ 
are of the same order 1. and s r k collect all the next-to- leading contributions, which are 
counted as 0(1/L 2 ). Consequently, for a fixed value of the condensate E, the inverse mo- 
ments (((o^) 2 ))^^ /V 2 and ((ui))^^ /V 2 can be expressed in the form Ylt=~i a nL 2n . The 
even powers n = 2,4 are the original leading terms, whereas the odd powers n = —1, 1,3 
arise from the next-to-leading corrections due to the non-zero modes. Hence, this type of 
correction does not mix with the leading contribution as far as the volume dependence is 
concerned. 

This is not true for s^, which stems from the zero- mode contribution of £3. They modify 
the constant term n = of the sum rules, and may be considered as small to the extent 
that E is small (let us recall that the dimensional estimate of the low-energy constants TZi 
leads to TZi ~ F 2 /A H with A H ~ 1 GeV). Of course, close to the critical point, one precisely 
expects E to become small. 

In the case of the A~ 6 -sum rules ( |4.33j )- (|4.35| ), the situation is similar, but now, the 



constants TZi already affect the coefficient of C 1 in the sum rules (they also change the 
constant term 
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VI. EXTRACTION OF PARTICULAR LOW-ENERGY CONSTANTS 



Near the critical point n cr it(N c ), we would like to exploit the sum rules for (((a 2 ) 2 )) ( Nf) 

and ({(Ti))^^ in order to isolate particular ratios of low-energy constants present in C 2 . in 
particular, it would be interesting to obtain a ratio with a specific sensitivity to the phase 
transition. To reach this goal, it is preferable to eliminate the next-to-leading corrections, 
which involve either unknown parameters like TZi or regularization-dependent quantities like 
g. As already pointed out, Eqs. (|5.26|) and (|5.27| ) can be viewed as expansions in the variables 
EL 2 and ML 2 . Hence, they are even functions of L. 



A. Varying the size of the box 

To exploit the structure of the sum rules at the next-to- leading order, it is therefore 
interesting to introduce the derivative-like operator: 



X 



{ (L - |) F(L + a) + (L + |) F(L - a) - 2LF{L)) 



where a is an arbitrary parameter. If we consider an even monomial F(L) = L 2k , S a [F] is 
an even polynomial of degree L 2k ~ A . We obtain for the first powers: 

L°^0, L 2 ^0, L A ^a\ (6.2) 
L 6 ^a 4 (3L 2 + 2a 2 ), L 8 -> a 4 (6L 4 + 12L 2 a 2 + 3a 2 ). (6.3) 

If we denote £ 2 = (((cx 2 ) 2 }} { Nf) /V 2 and t 4 = ((a A )) { " f) /V 2 , we have: 

S a [N f t 2 - U](L) = ^ 2z s + 0(1/L 6 ), (6.4) 

Sa[6aM(L) = 8Ar , ( 3 ^_ 1} S 4 + 0(1/L 10 ). (6.5) 

In order to get a quantity which is invariant with respect to the QCD renormalization 
group, we take the ratio of these two sum rules: 

The evaluation of p requires the knowledge of t 2 and £4 for five different box sizes: L — 2a, 
L — a, L, L + a, L + 2a. Notice that a is not required to be large; it is sufficient to have 
L — 2a much bigger than 1/A#. On the other hand, for too small a, the difference operator 
5 a may be too sensitive to numerical errors. 

For a discretized torus (a lattice), we can put L = na with n integer and a the lattice 
spacing. Eqs. (|6.4|) - (|6.6|) remain true, and the comparison of different volumes is translated 
into the evaluation of the inverse moments on lattices with the same spacing, but with 
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various numbers of sites. The powers in the lattice spacing a on the right-hand side of 
Eqs. (|6.4p -( |B~6|) reflect merely the dimension of the involved quantities. 

Eqs. (|6.4|)-(|6T6D are independent of the next-to-leading order contributions: this allows 
to consider smaller volumes than previously stated (for instance, the estimate L > 1.9/ F 
of Sec. |V A| , based on our next-to-leading order analysis does not necessarily apply to the 
sum rule ( |6.6|) ). The volume-independence of Eqs. (|6.4j )-( |BT6l) could already be seen for 
smaller volumes. Moreover, the inverse moments must fulfill another non-trivial consistency 
relation: 

p'= 6 j§^(L) = ± + 0(l/L™). (6.7) 

The ratio p is invariant under the QCD renormalization group and its variations with Nf 
could reflect the proximity of the critical point in a particularly clean way, as discussed in 
the next section. 



B. Relevance of the ratio 

We have argued in a previous paper || that the approach to n CI it(N c ) could result into 
a large Zweig-rule violation in the scalar channel. Let us recall briefly the argument. We 
consider the chiral limit for the first Nf light flavours of common mass m — > 0, and denote 
by s the {Nf + l)-th quark, whose mass m s is non-zero, but still considered as light compared 
to the scale of the theory (real QCD corresponds to Nf = 2). Here, we typically consider 
Nf such that Nf + 1 < n cri t(iV c ) < Nf + 2. E(iVj) is a function of m s , with the derivative: 

-^—E(Nf) = lim f dx(uu(x)ss(0)) c = U z (m 8 ), (6.8) 

om s m->0 J 

where the superscript c stands for the connected part. Since E(iVy) — > E(iV/ + l) for m s — > 0, 
one can write: 

Z(N f ) = E(N f + l)+ dp:n z (fi) =£(JV) + 1) +m s Z eS (m s ) + 0(m 2 s log m s ), (6.9) 

Jo 

where Z eff (m s ), defined in Ref. 0, is essentially 2Z s {Nf + 1), up to corrections of the 
order (S(A^j + l)) 2 , which are small in the vicinity of the critical point. Close to ^^(./Vc), 
the condensate term need not dominate the expansion (|6.9|) in powers of m s , due to the 
suppression of £(iV/ + 1). The large variation of the quark condensate from Nf to Nf + 1 
is then reflected by a large value of Zs{Nf + 1), related to the Zweig-rule violation in 
the ++ channel. Once expressed through the Dirac spectrum, E can be interpreted as the 
average density of small eigenvalues, whereas Hz is related to the density-density correlation. 
The ratio ZgjY? measures therefore the fluctuation of the quark condensate. For Nf near 
the critical point n crit (N c ) where E vanishes, one may expect a suppression of E and an 
enhancement of its fluctuations Z$- 

We can express the ratio Zs/T< 2 by introducing the Gell-Mann-Oakes-Renner ratio [|3T 
measuring the condensate in physical units: 
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XGOR[Nf) 



2m£(A / ) 
FZM1 ' 



(6.10) 



where m denotes the common mass of the Nf lightest quarks (m = (m u +md)/2 for Nf = 2). 
Following the analysis of Ref. ||, one obtains from Eq. ( |6.9p , in the approximation Z c g ~ 
2Z S : 



F*Z s (N f + 1) 
I?(jV) + 1) 



eff 



2E 2 (A / + 1) 

X GOR (iV / )-X G o R (iV / + l) 
[X GO R(N f + l)] 2 



2rM 2 



+ 



(6.11) 
(6.12) 



where r stands for m s /m and the dots denote higher-order terms. For Nf <C n crit (N c ), 
the right-hand side of Eq. ( |6.12|) is very small. It can be illustrated by choosing Nf = 2, 
X(2) ~ 0.9 and r ~ 26 (Standard estimates). The difference of the GOR ratios 

satisfies in this case the lower bound: X(2) — X(3) > 0.2 and we consider this 

bound conservatively as an equality. In this case, the right-hand side of Eq. (|6.12 ) is of the 
order of 10~ 2 (let us notice that in this case, this quantity is related to 16L 6 (/z) at a typical 
hadronic scale fi ~ M p , c.f. Ref. [^UJ). The proximity of a phase transition could be detected 
by a considerable increase of the ratio ( |6.11| ) compared to its typical size ~ 10 -2 . 



C. Application to the lattice 

An evaluation of the inverse moments through lattice simulation represents a few in- 
teresting features. We work at finite volumes: the volume dependence is crucial to obtain 
information on the relevant low-energy constants of the effective Lagrangian, and the ex- 
trapolation to an infinite volume is avoided. The limitation to the topologically trivial sector 
is natural on the lattice by choosing strictly periodic boundary conditions. 

We do not aim at solving full QCD on the lattice. We want to compute Dirac inverse 
moments, averaged over the gluonic configurations with the statistical weight ( |2.13|) . To 



perform this more limited task, we have to know the Dirac spectrum for each gluonic config- 
uration. It can be obtained through the square of the Dirac operator: Ip 2 = D 2 + iF^a^. 
It seems much simpler to discretize this operator instead of Tj) itself. In particular, the dou- 
bling problems are not expected to arise in the spectrum of an elliptic operator like D 2 . It 
should be stressed that, while this procedure could be applied in our particular problem, it 
can hardly represent a general solution for doubling in the spectrum of lattice fermions. 

For a given gluonic configuration, we can therefore compute the inverse moments from 
the Dirac spectrum (which is independent of the number of flavours). The essential con- 
tribution to each inverse moments stems from the lowest eigenvalues. In this case, the 
Ay-dependence in the average (())o is expected to arise mainly through the Ay-th power 
of the product of the lowest Dirac eigenvalues, i.e. from the infrared part of the truncated 
fermion determinant, c.f. Eq. ( 2.1 9|) . The ultraviolet part of the determinant should then be 



included by a matching with the perturbative tail as discussed in Ref. |32] . A first possibil- 
ity consists in generating the gluonic configurations in the quenched approximation, and to 
include explicitly the fermion determinant in the observable. The advantage of this method 
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is that it would allow to change easily and continuously Nf while keeping the same set of 
gluonic configurations. On the other hand, Monte-Carlo simulation of the pure gauge the- 
ory could lead to a rather different distribution of small Dirac eigenvalues than a simulation 
including the fermion determinant into the statistical weight: the quenched generation of 
the configurations may therefore lead to biased results, when we use these configurations 
to compute quantities including explicitly a fermion determinant as an observable. If this 
reweighting procedure turns out to be inefficient, the generation of the gluonic configurations 
would have to include the product of the lowest Dirac eigenvalues in the statistical weight. 
The configurations should be regenerated for each value of Nf. 

The computation of the ratio p seems particularly attractive on the lattice. We have to 
compare five different lattice sizes to calculate this ratio, invariant under the QCD renor- 
malization group and protected from next-to-leading order effects. When Nf increases, an 
enhancement of 1/p would clearly indicate the vicinity of the critical point n crit (N c ) where 
the condensate vanishes. 



VII. CONCLUSION 

Two descriptions of Euclidean QCD on a torus can be fruitfully matched : the first 
involves the spectrum of the Dirac operator whereas the second relies on the effective theory 
of Goldstone bosons. The spontaneous breakdown of chiral symmetry can be related to the 
large- volume behaviour of inverse moments of the Dirac eigenvalues, En>oV^«i averaged 
over topological sets of gluonic configurations. Because of their sensitivity to Nf, these 
inverse moments can be used to detect chiral phase transitions occurring when the number 
of massless flavours increases. 

The quark condensate T,(Nf) is the chiral order parameter that is the most sensitive 
to Nf. It is conceivable that just above the first critical point n cr i t (N c ) where £ vanishes, 
the chiral symmetry is only partially restored. Below this critical point, the large- volume 
behaviour of the inverse spectral moments is given by the Leutwyler-Smilga sum rules (|3.12|) - 
( p.!4j ) and it is driven by the quark condensate (this behaviour corresponds to eigenvalues 
accumulating as l/L 4 ). Above n crit (iV c ), the asymptotic volume dependence of the inverse 
moments changes (see Eqs. (|3.16| )- (|3.17| )), corresponding to eigenvalues behaving as 1/L 2 . In 



this case, the dominant contribution comes from terms in the effective Lagrangian quadratic 
in quark masses. 

When Nf increases and approaches n CTit (N c ), the quark condensate becomes small, and 
its fluctuations (related to the Zweig-rule violation in the scalar channel) are expected to 
become large: the terms of the effective Lagrangian linear and quadratic in the quark masses 
may therefore contribute with a comparable magnitude. Hence, it may become necessary to 
include both of them into the leading order of the expansion of £ e ff, in order to derive the 
large-volume behaviour of the inverse spectral moments, which interpolates between both 
phases. The resulting sum rules have been analyzed in the topologically trivial sector v = 
(see Eqs. ( 4.31 )-( 4\35|) ). In particular, the formulae concerning positive inverse moments 



restrict the parameters of the effective Lagrangian. 

For Nf <C n cr it(-/V c ), the first subleading corrections to Leutwyler-Smilga sum rules are 
due to the non-zero modes, and reduce to a volume-dependent redefinition of the low-energy 
constant E. The next-to- leading corrections to these formulae have been calculated also for 
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Nf close to the critical point n cr it{N c ). The part arising from non-zero modes is translated 
into a redefinition of the low-energy constants E, Z$, Zp and A. The NLO contribution due 
to zero modes can be computed directly for v = 0. All NLO corrections behave as 0{L~ 2 ) 
relatively to the leading contribution. 

We have shown that combining inverse spectral moments at different volumes allows 
one to isolate the ratio of low-energy constants i?s/E 2 which is particularly sensitive to the 
chiral phase transition. The resulting "five-volume formula" ( |6.6| ) is furthermore insensitive 
to NLO finite-size corrections, and it is invariant under the QCD renormalization group. 

The study of the inverse spectral moments of the Dirac operator seems a promising tool 
to investigate chiral phase transitions in association with lattice simulations. The sums 
over eigenvalues can be computed from a set of gluonic configurations with v = and the 
corresponding Dirac spectra, obtained after the diagonalization of D 2 + Fa/2. The Nf- 
dependence is explicit, via the infrared part of the truncated fermion determinant and the 
finite- volume effects are not only taken into account, but essential for our purposes. 

The possibility to vary on the lattice parameters fixed in the real world, like Nf (and 
N c ) could open a new window on the chiral structure of QCD vacuum. This investigation 
could lead to a better understanding of QCD-like theories in general. For instance, among 
electroweak symmetry breaking models, technicolor and similar proposals have often been 
ruled out, assuming a smooth and simple dependence on Nf and N c leading to a direct link 
with actual QCD phenomenology ||33|| . If the chiral phase structure of vector- like confining 
gauge theories turned out to be richer, the chiral symmetry could be broken following a 



different pattern from actual QCD, offering new possibilities for technicolor-like models p4 
The study of A^-induced chiral phase transitions could therefore represent a step towards 
alternative theories of electroweak symmetry breaking. 



ACKNOWLEDGMENTS 

We thank Ph. Boucaud, L. Girlanda, P. Hasenfratz, H. Leutwyler, G. Martinelli, 
B. Moussallam and C. Roiesnel for valuable discussions. 



APPENDIX A: INTEGRATION OVER UNITARY MATRICES 

In the flavour space, we can define a complete set {t a } of N 2 Hermitian Nf x Nf matrices 
generating U(iV/). a is an index from to N 2 — 1: to is proportional to the identity, and 
the other matrices are traceless. They are normalized by: 

(tah) = ~Sab, J2 tata= l N f> ( Al ) 

a 

with the interesting identities for any matrices A and B: 

J2(t a A)(t a B) = ±(AB), (A2) 

a 

J2t a At a = 1(A), J> a At aJ B) = \(A)(B). (A3) 
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We can decompose any complex matrix on this basis: X = ^ a X a t a . If we want to perform 
integrations over U(Nf) involving a unitary matrix U, the non-vanishing integrals have as 
many components from U = U a t a as from W = Y2 a U^a- The first ones are: 

[dU]= 1, (A4) 

U(iV» 

[dU]U a U; = ^-5 abl (A5) 

V(N f ) 

[ [dU]U a U;U c U* d = —L^(5 ab 5 cd + 6 ad 5 bc ) (A6) 

JV(Nf) iV / - 1 

16 

[du]u a u;ujjsujr f = 8 (a?) 

u(jv» N f{ N f -IJiNf -4) 

x {(^| - 2)[5 ab 5 c A/ + 8 ab 5 cf 5 ed + S ad 5 cb 6 ef 

+8ad$cf8 e b + ^af^cb^ed + ^a/^cd^efe] 
— 4 Ay [5 a fe (t c tdt e tf + t c tft e td) + S a d(tbt c tft e + t b t e tft c ) 

+8af{tbtctdte + tbt e tdt c ) + 5 cb (t a tdt e tf + t a tft e td) 
+5cd(tatbt e tf + t a tft e tb) + S c f(t a t b t e td + t a tdt e tb) 
-\-fieb{ta.tdtctf + t a tft c td) + 5ed{t a tbt c tf + t a tft c tb) 
+Sef{t a tbt c td + WdWft)] 
+ 16(t a tbt c tdt e t f + t a t b t e tdt c tf + t a tbt e tft c td + t a tbt c tft e td 

+t a tdt c tbt e tf + t a tdt e t b t c tf + t a tdt e tft c t b + t a tdt c tft e tb 

+t a tft c tdt e tb + t a tft e tdt c tb + t a tft e t b t c td + t a tft c tbt e td) 



APPENDIX B: LEADING-ORDER GENERALIZED LAGRANGIAN 

FOR TWO FLAVOURS 

For = 2, the situation is slightly different from the generic case, because SU(2) 
representations are pseudoreal. In particular, the correlator ((uu)[dd)), which defines Z s , 
contains a determinant-like invariant and is no more an order parameter. The leading order 
of the generalized Lagrangian for SU(2) is 

£ {2) = i {F 2 (2)(d^d li U} - 2£(2)(£/ t M + M*U) (Bl) 

-A(2){(U^M) 2 + (M^U) 2 ) - Z P (2){U ] M - M ] U) 2 
-H{2)(M ] M) -^ / (2)(detM + detM t )} . 

The new counterterm TC'{2) is consistently counted 0(p 2 ) in G%PT, since det M involves 
two powers of the mass. 
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Despite similarities between Ti'(2) and TC(2) (terms with no mesonic fields, absent from 
the low-energy processes), 7i'(2) is not necessarily divergent. In the Minkowskian metric, it 
can be defined through the chiral limit of the Zweig-suppressed correlator: 

2% J d 4 x e ipx (0\T{uu(x)dd(0)}\0) = H' r (2)(fi) + 0(p 2 ) Gx p T . (B2) 

It is easy to prove that, in the chiral limit, the identity operator, the quark condensate and 
the gluon condensate do not contribute to the Operator Product Expansion of ((uu)(dd)) Q, 
The correlator ( |B2| ) is dominated by d = 6 operators and it behaves as 0(l/p 4 ) for large 
momenta. It is therefore super convergent. 7i'(2) can be related to the scalar spectrum 
through a dispersion relation with no subtraction, similarly as in Ref. [ 2"0| . Despite the 
difficulty of estimating the resulting integral, 7i'{2) can be determined in principle from 
experimental data in the + sector, including not only the low-energy dynamics, but also 
information about higher resonances. 

Since 7i'(2) is free of ultraviolet divergences, we can formally rewrite the G%PT leading 
order of the two-flavour Lagrangian in the generic form (Nf > 3). We use the identity, true 
for any 2x2 matrix C: (C) 2 — (C 2 ) = 2detC. This leads to a formal identification: 

A(N f ) <- A(2) - H'{2)/2, Z s (N f ) <- W'(2)/4, (B3) 
Zp(N f ) <-> Z P {2) + W'(2)/4, H(N f ) <-> H(2), (B4) 

which enables us to treat the two-flavour Lagrangian in the same framework as the generic 
case, even though the phenomenological interpretation of its parameters is different. 

APPENDIX C: EXPANSION COEFFICIENTS OF THE PARTITION FUNCTION 

This section is devoted to the calculation of the coefficients arising when the partition 
function is expanded in powers of X = ML 2 for Nf near (but under) the critical point 
n crit (N c ). The main lines of the computation are exposed in Sec. |1V B| , but its technical 
details and the results for an arbitrary winding number are presented here. The coefficients 
a u , /3 V . . . are defined in Eq. QL8p. 

1. Leading coefficient a v 

To compute a v (b, z, a), we begin with a u (b, 0, a), given by the leading order in x of the 
group integral: 

J° =!„(&, a; X • 1) = / [dU}(detUyexp[bx(U f ) + ax 2 {U ]2 )}. (CI) 



"Basically, the quark condensate cannot appear in OPE of ( |B2| ) without a mass term, vanishing 
in the chiral limit, whereas the discrete symmetry ul — > —ul rules out the identity operator and 
the gluon condensate. We thank B. Moussallam for this remark. 
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We can use Weyl's formula to transform the group integral into an integration over the 
eigenvalues of U: exp(z0 fc ) (& = !... Nf): 



[dU] 



V(N f ) 



n£n*i 



. fe=i 



(C2) 



with P = Ylk<i ( e ^ fc — e 1 ^ 1 ). P is a linear combination of exp(i ^2 nk4>k), with integers, 
antisymmetric under the exchange of two angles, so that for k ^ I, n/- and ni must be 
different. Their set forms one of the Nf\ permutations of (0, 1, 2 . . . Nf — 1), and P collects 
all of them, with a sign depending on the signature of the permutation. If the integrand is 
symmetric under the angle permutations, PP* can be rewritten [|21||: 



IPI 



N f \ e(^)exp 



fc=l 



with V(Nf) the set of the permutations over (1 . . . Nf) and e the signature. 
The group integral I" becomes: 

(N f \ AT/ 

n -t |p|2 n «* [««-* + - 2 ^ 2ifc ] ) • 
fe=i / fe=i 



(C3) 



(C4) 



When \P\ 2 is replaced by its symmetrized value ( |C3"D , the integrals over the angles become 
independent of each other: 



N f 



e ^>n 

r&V(N f ) k=l 



2tt 



exp [h 



Nf 



E 'Mil**' E 



a£V(Nf) 



k=l 



p k +2q k =s(k) 
Pk,Qk>0 



<xe 



(C5) 
(C6) 



with s(k) = k — a(k) + v. Obviously, if s(k) < for at least one k, the permutation does 
not contribute. But V(Nf) includes the identical permutation and v > 0: there is at least 
one contributing term in J", and all these contributions lead actually to the same leading 
power in x: 



N f 

n 

k=l 



X' 



,(k) 



X 



E k-a{k)+v _ x vN 



(C7) 



which is consistent with the factor (detX)^ in the expansion ( f4.8|) . We get therefore: 

a u (b, z = 0, a) = b vN f- 2m a m c m , 

m=0...uN f /2 



(Q 
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with the purely combinatorial coefficients: 



E <°) E 

aeV(N f ) {q k =l...s{k)/2} 
E Qk=m 



l[q k \(s(k)-2q k )\ 



Another way to describe c m is the generating polynomial: 



m=0...uN f /2 



X u -N+1 Xv-isf+2 -X"i/-JV+3 



with the polynomials in w: 



X 



E 



' -2q)Y 

q=0...j/2 H KJ H > 



x v 


X u+ \ 


X u+ 2 


■ ■ X u+N _i 




x v 


X u+ \ 


■ ■ X u+N _ 2 


X u -2 


X v -\ 


x u ■ 


■ ■ X u+ N-3 



x v 



(C9) 



(CIO) 



(Gil) 



Since the derivatives of I v with respect to b and z are not independent, Eq. ( |4.7[ ) yields 
the general expression of a v : 



a v (b,z,a)= b l a m z p ' ,, , — <■„,. 

l+2m+2p=vN ; 



l\p\ 



(C12) 



2. Subleading coefficients (3 U , 7„, <5j, 



We denote the various derivatives of a v : 



a, 



a,. 



d 2 a u da u 



(C13) 



db da u db 2 dz 

For b = b, a = a, z = z, for N flavours, and denoting K — N + \v\, the coefficients are 



1 



P = a— (y + b 2 ) + a'— —b (2Nz + 2a + Ny) + a"——y(Nz + a) 



K 



KN 



NK' 



(C14) 



7 = a 



_L_[^ + 2b 2 y + 2b 2 z + y 2 + 2z 2 + 2a 2 ] 



K(K 2 



-2a[b 2 + 2z\ + 



(K - N)(KN + 1) 
K(K 2 - l)(N 2 - 1) 



2a z 



+a' 



1 



K 2 - 1 
1 



b[b 2 y + 2b 2 z + 2y 2 + 6yz + 4z 
1 



K(K 2 - 1) 



2ab[y + 2z] + 



N(K 2 



2ab[b 2 + 2y + 2z) 



(C15) 
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KN(K 2 -1) K(K 2 — l)(iV 2 — 1) 

+a" ^ J J— [ [h 2 y 2 + 46V + 2b 2 z 2 + 5y 2 z + 4z 3 } 



NK(K 2 - 1) iV(K 2 - 1) 

1 r 9 9l i(JV + 1 9l9 

a[y 2 + 4z 2 ] + ——————— —2a V 



X(X2_ 1 ) J jVA:(jV 2 - 1)(K 2 - 1) 

W" \j^—2byz[2z + y] + J ^—- ) 2aby[Az + y] 
KN + 1 2 

+ 7vx(£: 2 -i)(7v 2 -i) ^ 
WW {^ 2 ^ 22 + iv(^I) 4a ^ 

KN + 1 



NK(K 2 - 1)(N 2 - 1) 

[a2a 2 (2z + b 2 ) + d'4a 2 % + &"2a V] 



NK(K 2 - 1)(N 2 - 1) 
6 4 



a i"^C^T) [ 2" + 2}?z + 2b2y + y2 + 2z ' + 2 ° 2] 

1 r , 9 n (AT- N)(K + N) , 

+^T 2 ^ 2 + ^- V-i)(iv 2 -i) 2a 

+J^2ab[y + 2z\ - NK{ ^ 2 _ i) 2ab[2y + 2z + b 2 ] 
1 o 2, KN + 1 - 

+ iv(^-i) 2a 6 + ir(^-i)(iv 2 -i) 2a 6 

+a" {- ^^.^ t^V + 4 ^V + 2&V + 5y 2 z + Az 3 ] 

1 r 9 9n K + N 9 9 1 

+ ^T a[ ^ + ^ 1 - NK(N 2 -l)(K 2 -l) 2a b j 

W" {-^3Ty 2 ^[^ + y] - ^^y 2 ^ + y] 



NK(K 2 - 1)(N 2 - 1) 



4a ay 
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NK(K 2 — 1)(N 2 — 1) 

k A ' [a2a 2 (2z + b 2 ) + a'4a 2 by + a"2a 2 y 2 ] 



NK(K 2 - 1){N 2 - 1) 



APPENDIX D: DIMENSIONAL REGULARIZATION ON A TORUS 



Following the regularization procedure described by Hasenfratz and Leutwyler ||30|| , we 
want to regularize sums like: 



(Dl) 



where if is a function and p is summed over 2tt/L ■ Z 4 . The Fourier transform of H(p) is: 

d d p 



H x) 



(2tt)' 



and fulfills the identity: 



V 

V l 

where I is summed over L ■ Z 4 . Because of the relation: 
m Gr = lim — H(v) — I — 



(D2) 



(D3) 



it is possible to separate in Gh the cut-off and the volume dependences: 



Gh = lim G H + gH, 



(D4) 



(D5) 



The infinite- volume limit of Gh contains the divergences for d — *■ 4 and has to be regularized, 
whereas Qh depends only on the volume. 
For v = ^2 1/n 2 , we have the relations: 



1 1 



lim 



-y— 



V ^ p 2 + M 2 VM 2 





\r 1 1 


= lim 


H VM 2 _ 


Af-»0 



(D6) 



with H(p) = l/(p 2 + M 2 ). In the case of the dimensional regularization, (P5|) involves: 



M 2 1 8 

G H = —(lnM + Cl ) + g H , g H = ttt^ - t? + 0(M 2 



W L 2 



(D7) 



where Ci contains a pole for d = 4, and /?i is a constant called "shape coefficient" , depending 
on the geometry of the box. For a four-dimensional torus, (3i = 0.1405. The dimensional 
regularization yields finally: 
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v = V^--47r 2 A (D8) 

For u, we can follow the same guideline and take H{p) = 1 for p ^ and H(0) = 0. Its 
Fourier transform is H(l) = 5^(1). gn vanishes, and we know that dimensionally regularized 
integrals like f vanish as well, so that: 

m = ^1^0. (D9) 
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FIGURES 
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FIG. 1. Variations of R = ((c4))o/(((<72) 2 ))o as a function of the volume, measured in physical 
units F~ 4 (Nf = 3 flavours, Zp/A = —1/2 and ZgjA = 1/6). The variation of L§ is only a 
redefinition of the scaling parameter £ and leads to a global shift of the curves. The vanishing for 
intermediate volumes is commented in Sec. IV D| . 
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FIG. 2. Variations of R = ((o"4))o/(((o"2) 2 )}o as a function of the volume, measured in physical 
units F~ 4 {N f = 3 flavours, Z P /A= -1/2 and Z s /A = 1). 
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FIG. 3. Variations of S = (((ct 2 ) 3 )) /((ct 6 ))o clS EL function of the volume measured in physical 
units F~ 4 , for different values of S = 2s/ A (Nf = 3, L 8 = 0.1, P = —1/2). S is sensitive to the 
parameter S even for intermediate volumes. A different value of L% would merely lead to a global 
shift of the curves. 
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FIG. 4. Values of (Z S /A,Z P /A) for which the sum rule (H3l"l) for «(a 2 ) 2 ))J ) Ar/_3) 
for any positive scaling parameter £ = VY? / A (the forbidden zone is hatched). 
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FIG. 5. Values of (Z S /A,Z P /A), for which the sum rule (fO|) for ((a 4 )) 
any positive scaling parameter £ = VT, 2 /A (the forbidden zone is hatched) 
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FIG. 6. Values of (Z S /A,Z P /A), for which the sum rule ( gjgg ) for ((<r 6 )) 
any positive scaling parameter £ = T/£ 2 /.4 (the forbidden zone is hatched). 
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FIG. 7. Values of (Zs/A,Zp/A) for which the sum rule ( 4.32[ ) for ((o^))^ is positive for 



any positive scaling parameter £ = VT, j A (the forbidden zone is hatched). 
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